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We analyze the consequences that the choice of the output of the system has in the efficiency of 
signal detection. It is shown that the signal and the signal-to-noise ratio (SNR), used to characterize 
the phenomenon of stochastic resonance, strongly depend on the form of the output. In particular, 
the SNR may be enhanced for an adequate output. 



PACS numbers: 05.40. +j 

The phenomenon of stochastic resonance (SR) (l] Q 
has emerged in the last years as one of the most exciting 
in the field of nonlinear stochastic systems. Its impor- 
tance as a mechanism for signal detection has given rise 
to a great number of applications in different fields, as for 
example electronic devices [fjlj , lasers , neurons [ p^|Jl4| 
and magnetic particles [p.5|l6| . 

The most common characterization of SR consists of 
the appearance of a maximum in the output signal-to- 
noise ratio (SNR) at non-zero noise level, although dif- 
ferent definitions have been used in the literature. The 
definition through the SNR accounts for practical appli- 
cations, because the SNR is the quantity that gives the 
amount of information that can be transfered through a 
medium as well as measuring the quality of a signal. Ad- 
ditionally, the SNR quantifies the possibility to detect a 
signal embedded in a noisy environment. Another defi- 
nition of SR, apparently similar to the one of the SNR, 
has been proposed in terms of a maximum in the output 
signal. Although both the SNR and the output signal 
have been analyzed in terms of the parameters of the 
system, e.g. the frequency or the amplitude of the input 
signal, there is an important aspect which has not been 
considered in depth up to now. An adequate election of 
the output of the system may have implications in the 
behavior of the quantities used to manifest the presence 
of SR. This is precisely the problem we address in this 
paper. 

It is interesting to realize that normally the output of 
the system is the same as the dynamic variable x(t) enter- 
ing the stochastic differential equation, although some- 
times the sign function of x(t) has also been considered. 
No matter the system, the output may in general be any 
function of x(t) which is usually fixed through the char- 
acteristics of the problem. However, in order to detect 
a signal embedded in a noisy environment any function 
may be used. Thus, instead of Fourier transforming x(t) 
we can transform the function v(x(t)). 

Let us discuss one of the most simplest cases, in which 
the dynamics is described by an Ornstein-Uhlenbeck pro- 
cess, where the input signal modulates the strength of the 
potential in the following way: 



dx 
~db 



Here h(t) = k(l + asxn.(uiat)), with k and a constants 
and £(t) is Gaussian white noise with zero mean and sec- 
ond moment + r)) = D6(t), defining the noise 
level D. The effect of this force may be analyzed by the 
averaged power spectrum 
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dt 



(v{t)v(t + T))e-^ T dT . (2) 



To this end we will assume that it consists of a delta 
function centered at the frequency u>q plus a function 
Q(u>) which is smooth in the neighborhood of loq and is 
given by 



P(w) = Q[u) + S(u )6(uj - lj ) 



(3) 



Let us now assume the explicit form for the output of 
the system, v(x) — where j3 is a constant. Although 
this model does not exhibit SR, it is adequate to illustrate 
the form in which signal and SNR vary as a function of 
the output. Considerations about our model based upon 
dimensional analysis enable us to rewrite the averaged 
power spectrum as 

1 f D\ 

P(iv, D, k, a, uj , (3) = - I — J q{uj/uj , k/uj n ,a, (3) 

f)%(*Ma,^(l-^) , 

(4) 

where q(u>/ujo,k/u!o,a) and s(k/uo,a) are dimensionless 
functions. 

In spite of the simplicity of this result, a number of 
interesting consequences can be derived. From Eq. (Q) 
we can obtain the expression for the output signal 



S = 



s(k/wo,a,0) 



(5) 



— = -h(t)x+m 



(i) 



Three qualitative different situations are present depend- 
ing on the exponent 0. For /? > the signal diverges 
when the noise level D goes to infinity, whereas for (3 < 
the signal diverges when D goes to zero. Even more in- 
teresting is the case (3 — 0, in which the signal does not 
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depend on the noise level. The previous results have been 
verified numerically for some particular values of [3 (Fig. 
|l|), by integrating the corresponding Langevin equation 
following a standard second-order Runge-Kutta method 
for stochastic differential equations [ojlH]. 




FIG. 1. Output signal corresponding to Eq. ([l]) (k = 1, 
a — 0.5, and uo/2n = 0.1) for different exponents of the 
output (P = —0.5, 0, 1, 2, 5). The lines are fit to a power law 
Wfit = -0.501, 0.001, 0.998, 1.996, 4.993). 



SNR = k—— — — , (6) 

q(u;/uj ,k/uJo,a,(j) 

This result does not depend on the noise level thus in- 
dicating that the system is insensitive to the noise. No 
matter the noise intensity, the SNR has always the same 
value despite the fact that signal is a monotonic increas- 
ing or decreasing function of the noise. For further il- 
lustration of these features we have depicted in Fig. |^. 
the temporal evolution of the output of the system when 
v(x) = x 2 , for two values of the noise level. In both cases 
we have used the same realization of the noise. In the 
figure, we can see how the noise only affects the system 
by changing its characteristic scales. 

The former results refer to the behavior of the SNR 
as a function of D. For practical applications it is also 
interesting the knowledge of the SNR as a function of 
(3 based upon the possible increasing of the SNR when 
varying (3. We have found that the SNR has a maximum 
at j3 — 1 (see Fig. ||). Consequently, for the output class 
of functions v(x) = \x\° the input signal will be better 
detected when (3=1. 
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FIG. 2. Time evolution of x(t) 2 (same situation as in Fig. 
[j]) for the noise levels (a) D = 0.01 and (b) D = 1. 

It is interesting to point out that the signal increases 
for low or high noise intensities, depending on the value of 
the exponent [3. From the previous considerations it be- 
comes clear that the output signal itself does not always 
constitute a useful quantity to elucidate the optimum 
noise level to detect a signal. In contrast, the SNR over- 
comes this ambiguity. Its expression straightforwardly 
follows from Eq. (H) 




FIG. 3. SNR as a function of the exponent of the output 
(same situation as in Fig. |l|). 

All the functions we are considering as outputs are 
scale invariant, and dimensional analysis can be readily 
performed. However, when this requirement about v(x) 
does not holds, the previous results do not apply. This 
could be the case of the Heaviside step function v(x) = 
Q(x — 9), where 8 represents a threshold. In fact, this 
situation is quite similar to standard threshold devices 
H considered previously. In this case, both the SNR and 
the output signal exhibit a maximum at non-zero noise 
level (see Fig. ||). Although the evolution equation of the 
variable x(t) is linear, SR appears due to the fact that 
the output is a nonlinear function. 
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FIG. 4. (a) SNR and (b) output signal for the output 
equal to the step function with threshold — 1 (same situa- 
tion as in Fig. |lj). 

Having discussed the role played by the output in a 
simple monostable system, let us now analyze the case of 
the bistable quartic potential, which has been frequently 
proposed in order to describe the phenomenon of SR. The 
dynamics of the system is then given by the following 
equation: 

dx 

— = ax-bx 3 + A sin(w <) + (7) 

where a, b and A are constants and £(t) is the same noise 
as the one introduced through Eq. @). 

To study this system, one usually takes as output the 
variable x(t) and sometimes the sign function sgn(x(t)). 
In the limit when the amplitude of the input signal goes 
to zero, these two forms of the output give the same 
results (see ref. || for more details). However, when 
the input signal has a finite amplitude, the SNR for x(t) 
diverges, whereas for sgn(x(t)) goes to zero when the 
noise level decreases. Despite the divergence of the SNR 
for x(t), if the amplitude of the input signal is not large 
enough, the SNR has a maximum at non-zero D. 

As output, we could take in general x(t)P . The choice 
of f3 has important consequences as the SNR may depend 
on this parameter. Thus to better detect a signal, the 
noise level is not necessarily the only tunable parameter. 
In Fig. ||(a) we have plotted the SNR for different values 
of 0, observing its strong dependence on this parame- 
ter. In particular, for log 10 (-D) « —1.25, upon varying (3 
from 1 to 7 the SNR increases in about 12 dB. Moreover, 
when increasing (3 the maximum in the SNR becomes 



less pronounced and disappears for a sufficiently large (3, 
as occurs for the case (3 = 25 . In regards to the signal, 
variations of (3 change its behavior drastically. This point 
is illustrated in Fig. ||(b) where we can see that, when 
increasing the noise level, for (3 — 1 the signal goes to 
zero, whereas for the remaining cases the signal always 
increases at sufficiently high noise level. In Fig. |^(c) we 
have also displayed the output noise. From Fig. |^(a) 
it follows that a simple variation on the output changes 
the qualitative form of the SNR, in such a way that the 
maximum at non-zero noise level may disappear. Thus, 
the SNR is a monotonic decreasing function of D and ap- 
parently the input signal can always be better detected 
by decreasing the noise level. However, when the SNR is 
a decreasing function of D, there exists a region around 
the maximum, corresponding to the curves (3 = 1,3,5,7 
in which the SNR for (3 = 1, 3, 5, 7 is greater than that 
for (3 = 25. We then conclude that when increasing the 
noise level, the signal can be better detected if one simul- 
taneously changes the value of (3. 
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FIG. 5. (a) SNR, (b) output signal, and (c) output noise 
for the bistable quartic potential (a = 1, b = 1, ujo/2it = 0.1 
and A — 0.13) for the outputs f3 = 1 (empty circles), /3 = 3 
(empty triangles), /3 = 5 (empty squares), = 7 (filled cir- 
cles), and /3 — 25 (filled triangles). 

To end our analysis, in Fig. |^ we have displayed two 
temporal series for two different values of (3 at the noise 
level for which the effect of the variation on (3 is more 
pronounced. We can see that intrawell oscillations, for 
(3 = 7 are better observed than for /3 = 1. This fact 
explains the increase of the SNR. 
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FIG. 6. Time evolution of (a) x(t) and (b) x(t) 7 (same 
situation as in Fig. ||) for the noise level D = 0.056. The 
sinusoidal line in both figures indicates the value of h(t). 
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In summary, we have shown that the quantities (signal 
and SNR) used to characterize the phenomenon of SR 
strongly depend on the form of the output. In this regard, 
the behavior of the SNR has revealed to be more robust 
than the one corresponding to the signal. Our findings 
have important applied aspects since an adequate choice 
of the output of the system may be crucial in order to 
better detect a signal. 
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